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I. INTRODUCTION 



Good understanding of background and signal processes will be necessary to interpret 
data from the Large Hadron Collider (LHC) and observe physics beyond the Standard 
Model. In particular, large multiplicity final states are of interest Reliable predictions 
for such processes require computations of next-to-leading order (NLO) QCD corrections. 
Traditional methods for NLO calculations have difficulties in dealing with processes of such 
complexity; as a result, many new approaches to one- loop computations have been suggested 
in recent years 
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Among those approaches, genera 
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ized unitarity stands out 
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241 ] . The key feature of this method is that it 



allows calculation of one-loop scattering amplitudes directly from tree amplitudes leading 
to a computational algorithm of polynomial complexity 25[. The efficiency of generalized 
unitarity for NLO calculations for pro cesses with high multiplicity final states has been 
explicitly demonstrated in Refs. 26l. l27|. 

Until recently, generalized unitarity was mostly used to compute the cut-constructible 
of scatterin g am plitudes, while calculations of the rational parts proved to be 



parts 
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challenging. In Refs. 10, [ll| the four-dimensional boot-strap method was developed to 



evaluate the rational part. Another approach developed to generate the rational part uses 



generalized D-dimensional unitarity 



HQ. 



In a recent paper [24|, we extended the method of Refs. [23j, |29| in such a way that both 
cut-constructible and rational parts are obtained within a single formalism using integer- 
dimensional on-shell cuts. This method leads to a computational algorithm of polynomial 
complexity, as shown in Ref. 271 ]. 



Up to now, generalized unitarity has been mainly studied in the context of multi-gluon 
scattering amplitudes which simplifies the problem significantly. In the general case, one 
has to deal with two additional issues - different types of particles that participate in the 
scattering process and the fact that massive particles can be involved. It is necessary to 
address these issues before generalized unitarity becomes a practical tool for NLO calcula- 
tions of phenomenological interest. The goal of this paper is to do exactly that and extend 
the applicability of generalized D-dimensional unitarity by considering one-loop am plit udes 



involving gluons and massive quarks. The computational method developed in Ref. 24j can 
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handle both extensions easily. 

Dealing with particles of different flavors requires more sophisticated bookkeeping, but 
is otherwise straightforward. However, the presence of massive particles introduces new 
conceptual issues. An obvious consequence of having virtual particles with non-zero masses 
contributing to one-loop scattering amplitudes is that in addition to quadruple, triple and 
double cuts, we also have to deal with single-particle cuts. Such an extension is straight- 



forward; the necessary details have already been given in Ref. [23j. A more interesting 
consequence of massive particles present in the scattering process is that generalized unitar- 
ity applied to certain double- and single-particle cuts becomes more subtle. This is closely 
related to external wave function renormalization constants which originate from Feynman 
diagrams with self-energy insertions on external lines 1 . We will show that this complication 
can be circumvented without encumbering the formalism. 

To validate the method, we focus on the calculation of one-loop amplitudes with a mas- 
sive quark anti-quark pair and up to three gluons. These one-loop amplitudes have been 
calculated using more traditional methods. The one-loop corrections to tt + 2 and ti+ 3 
partons scattering have been first calculated in Ref. 30, |3l| and Ref. 32J , respectively. 

The outline of the paper is as follows. In Section II we discuss the modification of the 
D- dimensional generalized unitarity method required to include massive fermions. Section 
III describes the subtleties that arise when massive particles are involved in the one-loop 
scattering amplitude. In Section IV we present numerical results for the one-loop amplitudes 
t— > tt + 2 gluons and i — ^ -h 3 gluons. The conclusions and outlook are given in Section 
V. 



II. ONE-LOOP AMPLITUDES AND DIMENSIONALITY OF SPACE-TIME 

One-loop calculations in quantum field theory are divergent and require regularization at 
intermediate stages of the calculations. The conventional choice is dimensional regulariza- 
tion where momenta and polarization vectors of unobserved virtual particles are continued 



to arbitrary dimensions 
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34(. By keeping the momenta and polarization vectors of all 



observable external particles in four dimensions, one can define the one-loop helicity ampli- 



1 Similar problems appear due to diagrams that can be interpreted as one-loop expectation values of quan- 
tum fields. 
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tudes to be used in NLO parton-level generators 35|]. Once the dependence of a one-loop 
amplitude on the dimensionality of space-time is established, the dimensionality D can be 
interpolated to the non-integer value D = 4 — 2e. The divergences of one-loop amplitudes 
are regularized by the parameter e. 

While the analytical implementation of the dimensional regularization procedure is well- 
established (see for example Ref. 36j), a numerical implementation needs more considera- 



tion. In Ref. 24J we developed numerical implementation of dimensional regularization. To 
explain the method, we note that any iV-particle one-loop scattering amplitude A^j can be 
written as a linear combination of the so-called master integrals. The coefficients of such 
an expansion depend on D; this dependence can be made explicit by choosing the appro- 
priate basis of master integrals. After dimensional continuation, the final expression in the 
four- dimensional helicity (FDH) scheme 
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is given by 
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where we introduced the short-hand notation [ii|i r 
master integrals in Eq. (pQ) are defined as 

d D l 



AD) 

il—ir, 



1 < i\ < H < ■ ■ ■ < in < N. The 
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Tfi^ . The coefficients hi^i^^ bi^i^-i ^1*2*35 ^1*2*3 > ^*i?2*3*4' 



with dj = di(Z) = (Z +pi H pi) 2 

di^hn, di^izUi an d e hi2hui$ are independent of the dimensionality. 

We can compute these dimension-independent coefficients numerically, within the method 



of -D-dimensional generalized unitarity. To accomplis 
method 
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1 this, a parametric integration 



24j, based on the ideas developed in Ref. 



291 ] . is employed. The key point 



is to extend the dimensionality of the loop momentum to an integer -D-dimensional value. 
For one-loop calculations, an extension to five dimensions is sufficient 
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331 ] . However, 



care has to be taken with the dimensional dependence of the spins of the internal particles. 
The dimensional regularization scheme allows us to choose the dimensionality for internal 



degrees of freedom of virtual particles D s to be equal or larger than the embedded loop- 
momentum dimensionality. By choosing the parametric form of the integrand in integer 
(D s , D) dimensions we can determine the dimension-independent coefficients through par- 
tial fractioning of the integrand. The partial fractioning factorizes the calculation of the 



one-loop amplitude into tree amplitudes 24j. This factorization property is equivalent to 



the factorization obtained in generalized unitarity methods. The four-dimensional helicity 
scheme defines the parametric continuation as D s — > 4, D — > 4 — 2e with the constraint 
D s > D, giving the final result of Eq.flT]). 

We consider now the one-loop scattering amplitude involving a massive quark pair in 
addition to the gluons: i— > tt + N gluons. The -D s -dependence of the amplitude is linear. 
This means that we need to compute the integrand for two different values of D s so that we 
can separate the -D s -dependent and _D s -independent parts [24j. Because we need well-defined 



states for fermions when taking the internal fermion propagator on-shell, we must choose 
the space-time dimensionality to be even, i.e. D s = 4, D s = 6 and D s = 8. 

The on-shell internal gluonic polarization states in six and eight dimensions with the 
momentum vector in five dimensions are straightforward generalizations of the choices made 
in Refs. 0, 0]. The construction of /^-dimensional on-shell fermionic lines requires an 
extension of the four-dimensional Clifford algebra. We need to explicitly construct the D s 
Dirac matrices V 1 and the 2 Ds ^ 2 ~ 1 spin polarization states Uj(l,m) that satisfy the Dirac 
equation 

2 D s /2 D _ 1 

(E W - m x s *) u j\ 1 > m ) = °> ( 3 ) 

and the completeness relation 

2 (O s /2-l) D _ x 

where the on-shell condition for a fermion with the mass m and momentum I reads I 2 = m 2 . 
To construct the explicit higher dimensional Dirac matrices we follow the recursive definition 



given in Ref . 36| . The 8x8 six-dimensional Dirac matrices are defined in terms of the 4x4 



four- dimensional Dirac matrices {7 , 7 1 , 7 2 , 7 3 , 7 5 } 




-ii=l,2,3 
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It is readily checked that these matrices satisfy the standard anti-commutation relation 



ft" + ft" = 2g^, v = 0, 



(6) 



The 16 x 16 eight-dimensional Dirac matrices are constructed in a similar manner from 
the six-dimensional Dirac matrices. The D s -dimensional Dirac matrices are given for a 
particular representation of the Dirac algebra. Other representations can be obtained by 
unitary transformations. To construct a set of 2 Ds ^ 2 ~ 1 spinors satisfying the Dirac equation 
we generalize the procedures used in the four-dimensional case. We define the spinors 



u {s \l,m) 



For D, = 4 we choose 



v^o + m 
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1 2 D ^ 2 ~ l . 
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and construct recursively the D s = 6 eight-component basis spinors 
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The eight spinors for D s = 8 are obtained using the obvious generalization. It is easy to see 
that the spinors constructed in this way do indeed satisfy the Dirac equation. 

To check the completeness relation, we need the Dirac-conjugate spinor u. One subtlety 
associated with the fact that we have to deal with complex, rather than real, on-shell mo- 
menta is that in order to satisfy the completeness relation Eq. (j3J), we have to define the 
conjugate spinor as 

M(,, m)= ,M(i£+2>. (10 ) 

Note that the loop momentum is not complex conjugated. It is then straightforward to 
check that the completeness relation Eq.(j4j) is satisfied. 



III. MASSIVE PARTICLES AND THE UNITARITY CUTS 



To determine the dimension-independent master integral coefficients in Eq. (TjQ) we use 
the D- dimensional generalized unitarity method of Ref . 24J . To this end, we parameterize 
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the integrand of the one-loop amplitude 

r . ,. , di 1 di 2 di.,di.di 5 di-idi^di^di, ...... di,di 2 di„ _ di 1 di 2 d{, 

[n\*5\ [U|mJ lH |*3 J l»l|*2j Mm] 

(11) 

The left hand side of the equation is completely specified by the Feynman rules. The 
parametric form on the right hand side of the equation depends on a set of coefficients. To 
determine the coefficients for a given phase space point we use partial fractioning. This 
isolates the individual pole structures, thereby dividing the sets of linear equations to be 
solved into smaller subsets. More importantly, the partial fractioning sets groups of internal 
lines on-shell. This organizes the left hand side of the equation into products of gauge 
invariant tree amplitudes, thereby removing the necessity to compute individual Feynman 
diagrams to evaluate A$ (I) for a given loop momentum. 

This procedure can readily be applied in a situation when massive particles are involved 
in the scattering process. The presence of massive particles creates more types of master 
integrals or, equivalently, more different denominator structures in Eq. (Ilip . Furthermore, 
the single-cut (or tadpole) contributions to one-loop amplitudes have to be calculated so 
that the tadpole coefficient in Eq. (DD) can be determined. These issues complicate the 
bookkeeping, but do not add conceptual difficulties. 

However, a new conceptual issue does appear when dealing with the double cuts shown 
in Fig. [TJ Note that such cuts need only be considered for external massive states, since, if 
the external on-shell line carries a light-like momentum, the cut in Fig. [1] is set to zero in 
dimensional regularization. For massive particles these cuts do give non-vanishing contri- 
butions. The subtlety arising when such cuts are considered is related to a conflict between 
generalized unitarity and self-energy insertions on the external lines 2 . 

To see this, we study the contribution from a particular two-particle cut shown in Fig. [U 
The only outgoing external line to the left of the cut is the top quark and all other external 
particles are to the right of the cut. The residue of the one- loop amplitude for such a cut 
can be schematically written as 

Re S [A [1] (t, 9l ,...,g n ,t)]~ ^ A [0] (t,g*,t*) x A [0] (t*,g*,g h . . . ,g n ,t) , (12) 

states 

2 While we discuss the one-loop amplitude tt + N gluons, other processes with massive external lines can 
be treated in the same way. 
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FIG. 1: The external self-energy cut of a general ti + N gluon loop amplitude splits in the external 
self-energy contribution and the remaining higher point contributions. The different shadings of 
the blobs represent different content. 

where t* and g* denote the top quark and gluon cut lines respectively and the sum is over the 
intermediate states of the on-shell top quark and gluon particles of the two cut lines. The 
factorized on-shell tree amplitudes are given by A^(t,g*,t*) and A^(t*, g*, gi, . . . , g n , t). 
However, the latter amplitude is not defined. Separating the cut self-energy contribution as 
indicated in Fig. [1] gives for the tree amplitude 

M (t* , g*,gi,..., gn, t) = RP>9 m ,9i,~-,9n,i) + B{f% ^ ^ (13) 

{pt* +p 9 *r - mi 

Momentum conservation forces the invariant mass of t* + g* to be equal to the top quark 
mass squared, {p g * + p t *) 2 = ml making the one-quark reducible part of the amplitude 
singular. 

The singular contribution corresponds to the self-energy correction to the external top 
quark line. When one- loop scattering amplitudes are calculated using conventional Feyn- 
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man diagrams, these type of one-particle reducible diagrams are discarded; their effects on 
the scattering process are accommodated later through the external particle wave function 
renormalization constants. We would like to follow this approach in conjunction with the 
generalized unitarity technique, but then care has to be taken with the gauge invariance. 

Suppose we subtract the first term in Eq. fll3p from the tree amplitude; in recursive 
calculations this can be done by truncating the recursive steps. It is then easy to see that the 
remaining part of the amplitude B, the second term in Eq. ffl3l) . is no longer gauge invariant. 
Indeed, the discarded part of the amplitude is related to the self-energy correction on the 
external top quark line; such self-energy corrections produce on-shell mass and wave-function 
renormalization factors. While the mass renormalization constant, Z m , is independent of 
the gauge- fixing parameter, the on-shell wave-function renormalization factor, Z 2 , is not. 
For this reason we have to ensure that the gauge used in calculating the second term in 
Eq. (fT3|) and the gauge used in the calculation of the wave-function renormalization factor 
Z2 are the same. Since the wave-function renormalization factors are most easily computed 
in the Feynman gauge, we use this gauge to calculate the residue in Eq. (112j) . This means 
that the sum over gluon particle states for the cut in Fig. CD includes non-physical states 
such that 



Note that since the offending cuts never involve gluon self-couplings, ghosts do not need to 
be considered. 

Finally, we note that, for the most part, the coefficients in Eq. ( TTTT) are computed using 
the standard sums over physical states of the on-shell particles associated with cut lines. 
However, for a limited set of pole terms which contain the external self-energy contributions 
we need the procedure described in this Section. We emphasize that the conflict between 
unitarity and self-energy corrections to external massive lines is generic; it appears in any 
calculation of one-loop scattering amplitudes provided that massive internal or external 
particles are present. 




(14) 
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FIG. 2: The quadruple cut of the tt + 2 gluon amplitude decomposes into 3 gauge invariant 
contributions, each with its own 4-point master integral. The first box integral contributes to the 
primitive amplitude ^(lj, 2 i; 3, 4), the second to Al(1j, 3, 2 t , 4) and the third to 3, 4, 2 t ). 

IV. SCATTERING AMPLITUDES AT ONE-LOOP 



To implement the generalized unitarity method in a numerical algorithm, we decompose 



the i — ^ + 2,3 gluon amplitude into so-called primitive amplitudes |39| 3 . Within the 
context of D-dimensional unitarity the primitive amplitudes play a special role. Each prim- 
itive amplitude has unique unitarity cuts, i.e. the flavor of the cut lines is uniquely defined. 
This is shown in Fig. [2] for the example of the quadruple cuts applied to the i— > tt + 2 
gluon amplitude. This quadruple cut decomposes into three distinct gauge invariant cuts, 
each with its own master integral. Each of the three individual cuts contributes to one of 
the three primitive amplitudes A L (lf, 2 t , g u g 2 ), A L (l h gi, 2 t , g 2 ) and A L (l h g x , g 2 , 2 t ). 

The method described in this paper is amenable to straightforward numerical implemen- 
tation. To evaluate a primitive amplitude we consider all pole terms in the partial fractioning 
of Eq. ({TT]) . Double pole terms that correspond to massless two-point functions for light-like 
incoming momenta and single pole terms that correspond to massless tadpoles are discarded 



We adopt the conventions and normalizations of Ref. [39j to define the primitive amplitudes. 
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since the corresponding master integrals vanish in dimensional regularization. The tree am- 



plitudes for each cut are computed using Berends-Giele recurrence relations [40|. Because 
single particle cuts contribute, we need to evaluate the high multiplicity tree amplitudes 
it + it + 2, 3 gluons 

Before discussing numerical results for one- loop i— > ti + 2, 3 gluon amplitudes, we remind 
the reader that, when massive particles are involved, additional renormalization constants 
are required to arrive at physical predictions. In particular, for massive quarks, on-shell 
mass and wave function renormalization constants are necessary 4 . For consistency, we need 
those constants in FDH scheme. As described above, the wave function renormalization 
constant needs to be computed in the Feynman gauge. The bare quark mass m and the 
bare quark field ipo are renormalized multiplicatively 



m = Z m m, ip = V Z 2 ip. (15) 



We find (D = 4 - 2e) 



Z m = Z 2 = l-C F g s c r [ — \ l^^ + ^—)+°(9 s ,t) 



1 - CV^r ( 4V (- + 5)+0(gt,e) , 



where g s is the bare strong coupling constant, cr is the normalization factor, 

_ r(i + e )r(i - ef 

r (4vr) 2 -T(l - 2e) ' 1 ' 

C F = (N? - 1)/{2N C ) is the color factor and pL is the scale which is introduced in dimensional 
regularization to maintain proper dimensionality of the coupling constant. Note that in the 
last step in Eq. ( JT6l) we used D s = 4, as is required in the FDH scheme. 

We now present the results of the numerical evaluation of one-loop h- > ti + 2, 3 gluon 
scattering amplitudes in QCD. We do not include diagrams with closed fermion loops. In 
addition, external wave function renormalization constants and the coupling constant renor- 
malization factors are not included. However, we do include the mass counter-term diagrams 
which are necessary to obtain a result which is invariant under gauge transformations of the 



4 Note that the on-shell wave-function renormalization constant contains both ultraviolet and infrared 
divergences. Both show up as poles in e. 
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external gluons. For presentation purposes, it is convenient to normalize one-loop primitive 
amplitudes to tree-graph primitive amplitudes 



Af(lt,3,...,j -l,2 t , j,...,n) = c r I ^- + ^ + I A [ ?{l h 3, . . . , j - 1, 2 t , j, . . . , n 



a 



(i) u(j) 



A [ l\lt,3,...,n,2 t ) = c r 

jJCXl CX1J 

They can be extracted from Ref. 41| 



a 



L +^ + c ( ^|^ UJ (lf,3,...,n,2 t ). 



(18) 



The coefficients a^) and 6^ J parameterize divergences of the one-loop scattering amplitude. 
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The functions Sf it f i+1 = Sf i+1 j i depend on the flavor of particles /j , their momenta Pi 
and the scale \i. They read 
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Finally, we need to define the spin states of the gluons and top-quarks. For the gluons 
we use the conventional definition of the helicity vectors 



p^ = E(l, sin 9 cos (p, sin 6 sin </>, cos 6*) 

£ a{p) = cos ^ cos 4> T i si n 0; cos6> sin ± i cos (ft, — sin 6) . 

v2 

For the massive on-shell quarks (p = (E,p x ,p y ,p z ), p 2 = m 2 ) we use the spinors 

/ 1 \ / \ 



(23) 



u + (p) = vw+ 
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(24) 
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(25) 
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The numerical results reported below are obtained in conventional double precision us- 
ing a FORTRAN 77 program. The evaluation time does not depend on the helicities of 
the external particles but it does depend on the specific primitive amplitude. It takes less 
time to evaluate primitive amplitudes where quarks are adjacent, than to evaluate primitive 
amplitudes where quarks are separated by gluons. The reason for this is that it is computa- 
tionally more expensive to have more quarks involved in the evaluation of the primitive tree 
amplitudes. That is, the more quark propagators there are, the longer the evaluation time. 

For evaluating the master integrals we use the QCDLoop program developed in Ref. 42] . 
We have verified that our calculations correctly reproduce the divergent parts of primitive 
amplitudes, given in Eq. ([19]) . For all primitive amplitudes we have checked the gauge 
invariance by substituting a polarization vector of one of external gluons by its momen- 
tum. In addition we performed a Feynman diagram-by-diagram check on the results of the 
calculation. 



A. Scattering amplitudes with two quarks and two gluons 

In Table U we present the results for the finite parts of some of the primitive amplitudes 
A L ( If, 2 t ,g 3 ,gi), A L (lt,g 3 ,2 t ,g 4 ) and A L (lt, g 3 , g±, 2 t ) in the FDH scheme. The numerical 
results are obtained for the scale choice \i — E. We take the mass of the top quark to be 
m t = 1.75 and choose the following kinematic point (p = (E,p x ,p y ,p z )) 

p 1 = £(l,0,0,/3), p 2 = £ (1,0, 0,-/3), 

p 3 = £ (-1, sin 0,0, cos 0) , p 4 = £ (-1, -sin 0,0,- cos 0) • (26) 

with E = 10, — a/I — ml/ ' E 2 and = ir/3. Note that all the external momenta are taken 
to be outgoing. 

On a standard Pentium 2.33 GHz processor, it takes about 7.5, 11 and 12 ms to evalu- 
ate the primitive amplitudes Al(1j, 2 t , 3, 4), A^lf, 3, 2 t , 4) and A^lf, 3, 4, 2 t ) respectively. 
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Amplitude 


tree 




c 


+U +tj +3; +4 


0.0026595 i 


5.859738 +11.04762 i 


43.74436+11.04762 i 


+h +t, — 3, +4 


-0.127261 i 


18.49057-2.63910 i 


18.49058 -2.63910 i 


+h ~t, +3-, —4 


1.259555 i 


20.69972-0.144581 i 


20.52783 -0.14458 i 


+ti ~t, —3, +4 


-0.4198517 


22.16788 -3.40322 i 


22.68356 -3.40322 i 


+h +3, +U +4 


-0.0035643 i 


-0.26303343 


-0.26303305 


+U ~3, +t, +4 


0.170558 i 


15.2990066 


15.2990071 


+h +3; — U _ 4 


-1.688090 i 


20.8261462 


20.8261462 


+h — 3; — ti +4 


0.56269666 i 


22.0890527 


22.0890523 


+t) +3> +4, +t 


0.000905 i 


-26.24047 + 40.67377 i 


-123.4438 + 40.67377 i 


+ii ~3, +4, +t 


-0.043298 i 


20.00357 -1.69128 i 


20.00357 -1.69128 i 


+f> +3> _ 4 5 — t 


0.4285350 i 


21.83688-4.01097 i 


21.33165 -4.01097 i 


+?> — 3> +4, — t 


-0.142845 i 


19.48471 -1.85887 i 


21.00038 -1.85887 i 



TABLE I: The i— > it + 2 gluons results for the primitive tree amplitude and finite parts of 
the three one-loop primitive amplitudes for various helicities of gluons and top-quarks. Both the 
cut-constructible and total finite one-loop terms are given. 



Approximately half of that time is spent on the calculation of rational parts. We note that 

takes 



the calculation of the cut-constructible part of color-ordered four-gluon amplitude 
about 1 ms. The difference in CPU time between it + 2 gluons and the four gluon amplitude 
is not dramatic. The time difference is the result of several factors. First, a larger number 
of cuts has to be calculated. Second, in addition to the cut-constructible part we calculate 
also the rational part. Last, the evaluation of tree level amplitudes with (massive) quarks 
takes more computational effort. 



B. Scattering amplitudes with two quarks and three gluons 

In Table [TT] the results for the finite parts of the four primitive amplitudes 

A L (lt,2 t ,g 3 ,g 4 ,g 5 ), A L (lt, g 3 , 2 t , g 4 , g B ), A L {l h g 3 , g 4 , 2 t , g 5 ) and A L (lt, g 3 , # 4 , gs, 2 t ) are pre- 
sented in the FDH scheme. We take the mass of the top quark to be m t = 1.75, the scale 
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Amplitude 


tree 




c 


+h +t, +3, +4, +5 


-0.000533-0.000137 i 


9.584144+6.530925 i 


51.8809+6.543042 i 


+ii ~t, +3, -4, +5 


-0.004540 + 0.018665 i 


19.65913-11.77003 i 


23.00306-9.699584 i 


+ii +t, ~3, +4, ~5 


-0.004726+ 0.014201 i 


33.15950-1.832717 i 


33.71943 -3.142751 i 


+ii ~t, ~3, +4, +5 


0.045786 + 0.010661 i 


22.84043-6.540697 i 


23.03114-7.313041 i 


+h +3) +t, +4, +5 


0.000182 + 0.001369 i 


6.517366-1.277070 i 


19.37656+7.563101 i 


+h +3) ~t, —A, +5 


0.0467366-0.006020 i 


19.440997-7.639466 i 


20.93024-9.936409 i 


+ti _ 3, +t, +4, ~5 


0.019275 -0.0732138 i 


15.31910 -3.9278496 i 


15.176306-4.102803 i 


+ti _ 3, — i, +4, +5 


-0.018203-0.111312 i 


24.13158+1.431596 i 


24.70002+1.018096 i 


+i> +3) +4> +t> +5 


0.00060-0.001377 i 


13.13854+6.157043 i 


10.13113+13.83997 i 


+h +3) —4, — t, +5 


-0.047199-0.021516 i 


23.90539 -2.168867 i 


22.905695-4.284617 i 


+£> ~~ 3) +4> +t, ~5 


-0.015110+0.063118 i 


13.54258-7.800591 i 


13.50273-8.018127 i 


+ti ~ 3) +4> — t) +5 


-0.048800+ 0.112645 i 


21.77602+ 2.078051 i 


22.52784+1.424481 i 


+f> +3j +4> +5) +i 


-0.000252+0.000144 i 


-10.35085+45.26276 i 


-98.81384+52.81712 i 


+t> +3) — 4) +5) — i 


0.0050023+0.008871 i 


23.944473+2.862220 i 


20.92683-0.968026 i 


+h — 3j +4, —5, +f 


0.000561-0.004105 i 


-2.987822-42.00048 i 


-3.834451-43.67103 i 


+i! _ 3; +4) +5> _ f 


0.021216-0.011994 i 


19.72995-2.120128 i 


20.94996-1.684734 i 



TABLE II: The 1— > it + 3 gluons results for the primitive tree amplitude and finite parts of 
the four one-loop primitive amplitudes for various helicities of gluons and top-quarks. Both the 
cut-constructible and total finite one-loop terms are given. 

/1 = E and choose the kinematic point 

Px = E (1,0,0,/?), p 2 = E (1,0, 0,-/5), (27) 
p 3 = Eti (-1, 1, 0, 0) , p 4 = E£ (-V2, 0, 1, l) , 

P5 = ~Pl -P2—P3- Pi, 

where E = 10, /3 = ^l-m 2 t /E 2 and £ = 2/(1 + ^2 + ^) = 0.4823619098. 

On a standard Pentium 2.33 GHz processor it takes about 27, 35, 45 
and 50 ms respectively to evaluate the primitive amplitudes ^(lj, 2 4 , 3, 4, 5), 
A L (lt, 3, 2 t , 4, 5), A L (lt, 3, 4, 2 t , 5) and Al(1j, 3, 4, 5, 2 t ). Comparing these evaluation times 
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to the it + 2 gluon evaluation time, we see that the scaling is similar to the time scaling of the 
four and five gluon evaluation time in Ref. {23]. Similar to ti+ 2 gluons case, approximately 
half of the time is spent on the evaluation of the rational part. 



V. CONCLUSIONS 



In this paper we extended the method of generalized D-dimensional unitarity by com- 
puting one-loop scattering amplitudes for processes with massive quarks. We have proposed 
a solution to the subtleties associated with external self-energies and renormalization. We 
validated the method by computing the one-loop amplitudes for 1 — ^ tt + 2 gluons and 
t— > ti + 3 gluons. We have shown that the method is amenable to efficient numerical 
implementation. The results of this paper show that the generalized D-dimensional uni- 
tarity is a robust computational method. It will allow us to carry out NLO calculations 
for a large number of high multiplicity processes with massive particles, relevant for LHC 
phenomenology. 
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